GaN/AlN Quantum Dots for Single Qubit Emitters 
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We study theoretically the electronic properties of c-plane GaN/AlN quantum dots (QDs) with 
focus on their potential as sources of single polarized photons for future quantum communication sys- 
tems. Within the framework of eight-band k-p theory we calculate the optical interband transitions 
of the QDs and their polarization properties. We show that an anisotropy of the QD confinement 
potential in the basal plane (e.g. QD elongation or strain anisotropy) leads to a pronounced linear 
polarization of the ground state and excited state transitions. An externally applied uniaxial stress 
can be used to either induce a linear polarization of the ground-state transition for emission of single 
polarized photons or even to compensate the polarization induced by the structural elongation, 
(accepted at lOP Journal of Physics: Condensed Matter; ©lOP 2008) 

PACS numbers: 



I. INTRODUCTION 



Electrically triggered sources of single photons are es- 
sential building blocks for future quantum communica- 
tion and quantum computing systems. Devices based 
on single quantum dots (QDs) are promising candidates 
for reliable low-cost solutions in this matter. Single- 
photon emitters (SPEs) with extraordinary spectral pu- 
rity have been realized recently using InAs/GaAs QDs 
[l|, 0]. Here, emission with controlled linear polariza- 
tion is achieved at low temperatures by exploiting the 
fine-structure splitting (FSS) of the exciton ground-state. 
At elevated temperatures, however, this splitting, of 
up to 500 /xeV for large QDs Q, is smaller than the 
homogeneous line width, inhibiting a spectral separa- 
tion of the co-polarized line pair. Nitride-based QDs 
present important alternatives for single-photon emis- 
sion. In /i-photoluminescence (PL) and cathodolumines- 
cence (CL) spectroscopy narrow emission lines from sin- 
gle InGaN/GaN 0, i, g 0, [Eijlp, EI and GaN/AlN 
QDs grown on the c-plane [ij, [ij, [ij, [l^, [lg| and a- 
plane [l^ have been demonstrated. Also zinc-blende 
GaN/AlN QDs are under investigation [li,[ii|. In the fu- 
ture the spectral range of nitride QDs can be extended to 
the long-wavelength range for fibre-optical communica- 
tion employing InN/GaN [20|. Optically pumped single- 
photon emission from single GaN/AlN QDs has already 
been demonstrated [14]]. For InGaN/GaN QDs a polar- 
ization mechanism based on valence-band (VB) mixing 
effects has been discovered recently [111 . |2]| that leads to 
a separation of differently polarized lines more than an 
order of magnitude larger than the FSS in InAs/GaAs 
QDs. Similar polarization effects have been reported for 
c-plane GaN/AlN QDs shortly after ^^. 

In this paper, we show that VB mixing effects can be 
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exploited to control the polarization of photons emit- 
ted by single GaN/AlN QDs. First, the polarization 
mechanism observed for InGaN/GaN QDs will be re- 
visited in Sec. |TT1 Then, we will turn our attention to 
the GaN/AlN system. We calculate the polarization of 
the optical interband transitions for GaN/AlN QDs using 
strain-dependent eight-band k-p theory (Sec. IIIip . The 
model accounts for piezo- and pyroelectric effects and is 
implemented for arbitrarily shaped QDs. It is described 
in more detail in appendix \K\ 

We show that an asymmetric strain field in the growth 
plane inside the QDs leads to a linear polarization of 
the confined A- and B-type excitonic states in orthogo- 
nal directions. The in-plane strain anisotropy can orig- 
inate from different sources, such as QD elongation, in- 
homogeneous composition profiles, etc. As examples a 
slight elongation of the QDs (Sec. IIIIBp and, as an ex- 
situ approach, the application of external uniaxial stress 
(Sec. nil Cp are considered in detail. 



II. POLARIZATION EFFECTS IN INGAN/GAN 
QUANTUM DOTS 

We have shown recently that emission lines from single 
InGaN/GaN QDs are systematically linearly polarized 
in the orthogonal crystal directions [1120] and [TlOO] — 
a symmetry that is non-native to hexagonal crystals 
[111 . l2l| . The polarization of the emission lines has been 
attributed to the character of the hole state [formed by 
either the A or B valence band (VB)] involved in the re- 
combination process and a slight elongation of the QDs. 

The InxGai_xN/GaN samples were grown in the 
wurtzite phase by metal-organic chemical vapor depo- 
sition on Si(lll) substrates. The InGaN layers were 
grown at 800 °C with a nominal thickness of 2nm us- 
ing trimethylgallium, trimethylindium, and ammonia as 
precursors. The QDs are formed by alloy fiuctuations 
within the InGaN-layer [6]. The sample has been inves- 
tigated with single-QD CL employing Pt shadow masks 
to probe a few QDs only. Further discrimination of the 
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FIG. 1: (Color online) CL spectrum of a single InGaN/GaN 
QD. Grey lines belong to a different QD. The inset shows the 
evolution of the intensities of the five lives as a function of 
the polarization angle. 



single-QD lines was possible using their temporal jitter: 
All lines originating from the same QD show the same 
characteristic jitter pattern. Line groups of up to five 
narrow lines with the same jitter pattern could be iden- 
tified [ll|. An example of such a single-QD spectrum is 
shown in Fig. 1. All emission lines show a pronounced 
linear polarization. The polarization directions scatter 
around [TlOO] and [1120]. Both directions were found in 
each line group 11]. The emission patterns differ sig- 
nificantly from the well-known polarized emission spec- 
tra of II/VI and III/V QDs with cubic crystal structure. 
Exchange-interaction induced FSS of the exciton (see, 
e.g., i3, 23, 24]) or other excitonic complexes has not yet 
been reported for nitride QDs. 

In contrast to InAs/GaAs QDs or GaN/AlN QDs, 
which grow in the strain-induced Stranski-Krastanov or 
Vollmer- Weber growth modes, the InxGai_xN/GaN QDs 
form within the InxGai_xN layer due to composition fluc- 
tuations within the layer. For a theoretical description 
these QDs have been modeled as ellipsoids with linearly 
graded indium concentration [111 l2ll . |25| . The maximum 
indium concentration Xc at the QDs' centers is 50 %, the 
minimum concentration a^e at the QDs' borders is 5%. 
The QDs are embedded in an Ino.05Gao.95N quantum 
well with a height of 2nm. The QDs have a height of 
dz = 2nm and lateral extensions of d^/y = 4.6 — 5.8 nm. 
Henceforth x will denote the direction of the QDs long 
axis. Further details on the calculations can be found in 
appendix |^ 

The emission lines can be attributed to transitions in- 
volving either the confined hole ground state (ho = Ha), 
which is predominantly formed by the A band, or the 
first excited hole state (hi = hs), which is predomi- 
nantly formed by the B band. If a QD is slightly elon- 
gated, transitions involving Ha are polarized parallel to 
the elongation and transitions involving Kb perpendicu- 
lar to it [11|, [21|. The single-particle electron and hole 
orbitals are depicted in Fig. 2 for a circular QD (a) and 
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FIG. 2: (Color online) Confined single-particle electron (blue) 
and hole (yellow) states in a round (a) and an elongated (b) 
QD. The corresponding interband transitions are shown on 
the right-hand side. For the round QD, the A- and B-band 
transitions emit unpolarized light (with k \\ C), while the re- 
spective transitions in the elongated QD are polarized parallel 
and perpendicular to the long axis of the QD. The transitions 
between 3.33 and 3.36 eV originate from higher excited elec- 
tron and hole levels. 
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FIG. 3: (Color online) Degree of linear polarization of 
the eo-ho (black line and circles) and eo-hs (red line and 
squares) transitions as a function of the in-plane elongation 
of InxGai.xN/GaN QDs. 



a slightly elongated QD (b). The electron ground-state 
(eo) envelope functions have s-like symmetry; the enve- 
lope functions of the first excited electron state (ei) is 
p-like. In contrast, the first two hole states, Ha and hs, 
both have s-like envelope functions. Both states have siz- 
able oscillator strengths with the electron ground state, 
but behave differently if the strain field within the QD is 
anisotropic in the basal plane, e.g., due to an elongated 
QD shape: Ha aligns parallel to the QD's long axis, hs 
perpendicular to it. An analysis of the projections of 
both hole states on the k-p basis states Px and Py yields 
that the Px (Py) projection of Ha [hs) increases with 
increasing QD elongation, while the Py {Px) projection 
decreases [2l|. Consequently, the optical transition be- 
tween Ha {Jib) and eo is linearly polarized parallel (per- 
pendicular) to the long axis of the QD. The polarization 
increases with increasing QD elongation (Fig. [3]). 
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FIG. 4: Model structure of c-plane GaN/AlN QDs: Truncated 
hexagonal pyramid with height h between 1.5 nm and 3.5 nm. 
The lateral diameter d is determined by the dot's aspect ratio 
h:d = 1:10 (1:5). The wetting-layer thickness is w = 0.25 nm. 
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FIG. 5: (Color online) Radiative lifetimes of confined excitons 
in c-plane GaN/AlN QDs: Red lines and squares: Calculated, 
QDs with an aspect ratio of h:d — 1:10 . Blue line and circles: 
Calculated, QDs with an aspect ratio of h:d = 1:5. Black 
crosses: Experimental values from Bretagnon et aJ.[29l]. Black 
plus: Experimental value from Kako et al[14|. 



III. GAN/ALN QUANTUM DOTS 
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FIG. 6: (Color online) Energies of the single-particle electron 
and hole states. The electron states are labeled according to 
the symmetry of their envelope functions (s, p, and d shell). 
Hole-state labels are given according to the valence band from 
which the hole states are predominantly formed. 



trie fields within GaN/AlN QDs strongly affect the emis- 
sion energies and radiative lifetimes of localized excitons 
within the QDs via the quantum-confined Stark effect 
(QCSE) [13, In, H mi. In the center of the QD these 
fields are as large as 9.5MV/cm for QDs with an aspect 
ratio of 1:10 and 8.0 MV/cm for QDs with an aspect ratio 
of 1:5. Depending on the size of the QDs the radiative 
lifetimes range from a few nanoseconds for small QDs up 
to as long as 100 /iS for large QDs. Radiative lifetimes in 
this range would limit cut-off frequencies of devices to a 
maximum of a few kHz. Thus, with the fields present, 
only small c-plane QDs, i.e. those with high transition 
energies, are good candidates for fast single-photon emit- 
ters (SPEs). 



In contrast to the InxGai_xN/GaN QDs discussed 
in Sec. HIl GaN/AlN QDs grow strain-induced in the 
Stranski-Krastanov or VoUmer- Weber growth mode. 
Also their typical size and shape is known with greater 
accuracy. Experimental reports on the structural prop- 
erties of c-plane GaN QDs (e.g., Refs. [H [H, [H, [23, Hi) 
agree on the shape of the QDs, which is a truncated 
hexagonal pyramid. The reported heights (h) scatter be- 
tween 1.3 nm and 5nm. The aspect ratios {h:d, where 
d denotes the lateral diameter) in most reports are in 
the range of 1:5 to 1:10. The model structure derived 
from these reports is depicted in Fig. U) The model QDs 
have an aspect ratio of either h:d = 1:10 or 1:5. The 
height of the QDs is varied between 1.5 nm and 3.5 nm 
in steps of 0.5 nm. The thickness of the wetting layer is 
assumed to be w = 0.25 nm. This set of model QDs cov- 
ers the major part of experimentally reported QD struc- 
tures and yields excitonic transition energies and radia- 
tive lifetimes in good agreement with experimental values 
[ij, l2y| (Fig. O . The huge built-in piezo- and pyroelec- 



A. Single-particle Energy Levels 

The bound hole states in GaN/AlN QDs are — as in the 
case of lUxGai-xN/GaN QDs — formed predominantly by 
the A and B band. C-band contributions are small be- 
cause the biaxial strain within the QDs shifts this band 
to much lower energies [111, l25| . As a first approximation, 
for each band we expect to find a ground state with an 
s-shaped envelope function, which is only spin degener- 
ate. The p shell consists of two states and the d shell of 
three. Due to the different parities of the bulk conduc- 
tion and valence bands, the electron and hole states have 
a finite optical matrix element if their envelope functions 
have the same parity, i.e. the allowed transition chan- 
nels are s-s, p-p, s-d, etc. Each transition channel exists 
twice, once for the A-type holes and once for the B-type 
holes. Figure [S] shows the calculated single-particle elec- 
tron and hole energy levels of all QDs considered in this 
work including the six energetically lowest (highest) elec- 




unstrained 

--\B\)=ia-^^(\PA) + 'i\Pyi))-'ih\P..]) 
biaxially strained 

(e„. = £;» = -0-01) 

'~-'\A\)=^_i\PA) + AP„\)) 
asymmetricaliy strained 

(e,,,. = -n.01,e,„, = -0.02) 

'-~\A]) = \P,]) 
'- \P\) = \P,i^) 



FIG. 7: (Color online) Valence band structure of GaN in the 
vicinity of F. Upper panel: unstrained, center: under sym- 
metric biaxial strain in the basal plane, typical for c-plane 
QDs. Lower: With an additional strain anisotropy in the 
basal plane, e.g., in asymmetric QDs. The splitting between 
the two top bands increases and their characters change from 
A- and B-like to Pa;-like and Pj,-like. 



tron (hole) states. The electron states are all formed pre- 
dominantly by the conduction band (w 95 %). Therefore, 
the s, p, and d shells can be clearly distinguished. The 
hole spectra, however, are more complex, because holes 
states are not formed by either the A or B band, but 
by a mixture of both bands and even a small C-band 
contribution. Still we can characterize each hole state by 
the band that contributes the most to it (see labels in 
Fig. [H]). For QDs with an aspect ratio of 1:5 [Fig. EJa)] 
the A-band s-state (/iq = /i^; w 95 % A-band projection) 
and the B-band s-state {hi = /is; w 90% B-band pro- 
jection) are energetically well separated from the excited 
hole states. Both have an unambiguously s-shaped en- 
velope function (not shown here). The splitting between 
both states (w 9 — lOmeV) does not increase for smaller 
QDs, but is constant. It corresponds to the energy sep- 
aration between A and B band in strained GaN. The 
higher excited hole states can not be unambiguously as- 
signed to p- or d-like orbitals. Please note, that although 
they have been labeled according to the major band con- 
tributions, this contribution sometimes does not exceed 
50 %. The QDs with aspect ratio 1:10 [Fig.gljb)] exhibit a 
significantly smaller excited-states splitting for electrons 
and holes due to the weaker lateral confinement. An 
exception is the splitting between Ha and hB^ which is 
largely independent from the QD dimensions as discussed 
before. 



B. Optical Transitions and their Polarization 

In bulk GaN an anisotropy of the strain field in the 
basal plane changes the character of the band-edge states 
from A (B) type to P^ {Py), if the stress is larger along 
y (see Fig. [7]) . The optical matrix elements involving the 
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FIG. 8: (Color online) Oscillator strength between the single- 
particle electron and hole state for a QD with an aspect ratio 
of h:d =1:5 and a height of /i = 2.0 nm. Shown are all 
transitions involving one of the first six electron levels (ei) 
and one of the first twelve hole levels (hj), but only up to an 
energy difference et — hj of 3.75 eV. (a) For a symmetric QD. 
(b) For a QD with a 10% in-plane elongation along [1120]. 



confined A- or S-band hole states in the QDs respond 
differently to an asymmetry of the confinement poten- 
tial in the basal plane: They become linearly polarized 
in orthogonal directions. Figure[Ha) shows the oscillator 
strengths between the six lowest electron levels Ci and 
the twelve highest hole levels hj (up to an energy dif- 
ference Ci — hj of 3.75 eV) for the QD with h:d =1:5 
and h = 2.0 nm. Figure [H]Jb) shows the respective spec- 
trum for a similar QD with an in-plane elongation of 10 % 
along [1120]. Almost all transitions of the elongated dot 
show a pronounced linear polarization either parallel to 
the elongation or perpendicular to it. The transition be- 
tween the electron ground state eo and hA {hs) is lin- 
early polarized along (perpendicular to) the long axis 
of the QD. Both transitions together can be regarded 
as s-channel. The wave functions of the higher excited 
hole states are build of sizeable contributions from more 
than one valence band. Therefore not all of these states 
can be unambiguously assigned to p- or d-like symmetry. 
The p- and d-shell labels in Fig. [H] are given according 
to the main contribution. The QDs with aspect ratio 
h:d — 1:10 show the same polarization as those with 1:5 
albeit less pronounced. The strain anisotropy that is in- 
duced by the elongation affects the confined states less in 
QDs with larger diameter, because the main parts of the 
wave functions are located in the center of the dot. The 
energetic separation between the orthogonally polarized 
(coj/ia) and (eg, /is) transitions in the s-channel is inde- 
pendent of the size of the QDs or their vertical aspect 
ratio. For (in-plane) symmetric QDs it is about lOmeV 
and increases for anisotropic QDs. This large separa- 
tion between the orthogonally polarized lines should en- 
able a spectral separation of both lines even at elevated 
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FIG. 9: (Color online) Polarization degree of the excitonic 
ground-state transition as a function of uniaxial stress along x 
([1120]), relaxation of the sample in y and z direction accord- 
ing to the laws of continuum mechanics has been included. 
The calculations have been performed for three different QDs 
with an aspect ratio of h:d — 1:10 and a height oi h — 2.0 nm. 
Black Imes and circles: Symmetrical QD. Blue lines and di- 
amonds: QDs with 10 % elongation along x. Red line and 
squares: QDs with 10 % elongation along y. 



temperatures. Thus, the well-defined polarization of the 
ground-state transition can be exploited for polarization 
control in future SPEs. 



C. Control of the Polarization by Uniaxial Stress 



involving either A- or B-type hole states are polarized 
in orthogonal directions. The separation of the ^-type 
ground-state and the orthogonally polarized i3-type first 
excited state is ~ 10 meV and largely independent from 
the QD size and shape. 

A sufficient anisotropy can be induced either by a 
structural elongation of the QDs or by an externally ap- 
plied uniaxial stress. An in-plane elongation of the QDs 
of only 10% leads to a polarization degree of the exci- 
tonic ground-state transition of up to 6:1 depending on 
the other structural parameters of the QDs. An exter- 
nally applied uniaxial stress of about 300 MPa leads to 
a polarization degree of more than 5:1; larger stress re- 
sults in a complete polarization of the emission. More- 
over, a polarization resulting from structural elongation 
of the QDs can be compensated by external stress. This 
effect could be exploited for future devices, in particular 
to achieve a well-defined polarization in QD-based single- 
photon emitters. 
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Since the linear polarization of the transitions is a re- 
sult of the strain anisotropy, it can be induced by uniaxial 
stress, without any structural anisotropy of the dot. The 
effect of uniaxial stress on the polarization of the ground- 
state transition of the QD with aspect ratio 1:10 and 
height 2 nm is shown in Fig. [5] (black line and circles). A 
pronounced polarization is found at stress levels that can 
easily be induced by anisotropic strain relaxation in epi- 
taxial heterostructures (due to the formation of cracks or 
defects) or by externally applied stress. Thus, it is possi- 
ble to control the polarization of the emitted photons. 
Figure [HI also shows similar calculations for elongated 
QDs (blue lines and diamonds; red lines and squares). 
The polarization induced by the structure can be easily 
enhanced or even inverted by the external stress. 



IV. SUMMARY 

We have theoretically investigated the spectroscopic 
properties of c-plane GaN/AlN QDs with a special em- 
phasis on their suitability as sources of single polarized 
photons. We have shown that the linear polarization 
of the interband transitions is effectively controlled by 
an anisotropic strain field within the QDs. Transitions 



APPENDIX A: CALCULATION METHOD 

For the electronic structure calculations we use a three- 
dimensional eight-band k • p model implemented on a 
finite differences grid. The model accounts for strain 
effects, piezoelectric and pyroclectric polarization, spin- 
orbit and crystal-field splitting, and coupling between 
the valence bands (VBs) and the conduction band (CB). 
The k • p method is based on the one originally intro- 
duced by Stier et al. [32]. The implementation for 
wurtzite crystals is described in detail in Ref. 125|. De- 
tails of the calculation of radiative lifetimes are given 
in Ref. '33|. Few-particle corrections have been included 
using a self-consistent Hartree (mean-field) scheme j25| . 
For the InxGai_xN/GaN system material parameters 
from Ref. ^ have been included as described in Ref. ll|. 
Most of the material parameter used for the GaN/AlN 
system are based on the recommendations of Vurgaft- 
man and Meyer [S^]- For the parameters rnj,, Ai, and 
Ep we use recently developed band-dispersion parame- 
ters [3J, |3a] derived from accurate quasiparticle energy 
calculations [33] • For the valence-band offset between 
GaN and AIN we use the value of 0.8 eV determined by 
Wei and Zunger [3g] . As already detailed in Ref. [2^, we 
neglect the anisotropy of the static dielectric constants e,- 



and use the mean value for all direction. For the refrac- 
tive index of AIN we use the Sellnieier-formula derived 



by Antoine- Vincent et al[3£ 
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